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Abstract 

The phase ordering kinetics of the two-dimensional uniaxial nematic 
has been studied using a Cell Dynamic Scheme. The system after quench 
from T=oo was found to scale dynamically with an asymptotic growth law 
similar to that of two-dimensional 0(2) model (quenched from above the 
Kosterlitz - Thouless transition temperature), i.e. L(t)~(t/ln(t/t )) 1/ ^ 2 (with 
nonuniversal time scale t ). We obtained the true asymptotic limit of the 
growth law by performing our simulation for sufficiently long time. The 
presence of topologically stable 1/2-disclination points is reflected in the 
observed large-momentum dependence k~ 4 of the structure factor. The 
correlation function was also found to tally with the theoretical prediction 
of the correlation function for the two-dimensional 0(2) system. 
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The phase ordering of various systems with scalar, vector and tensor order param- 
eters has gained considerable interest over last few years |Tj. The system quenched 
from a high temperature homogeneous disordered phase into an ordered phase does 
not get ordered instantaneously, instead the various degenerate ground states com- 
pete to be selected [HE]- I n the process, the system develops length scales that grow 
with time and topological defects, if present, are eliminated. An infinite system will 
never achieve complete ordering and the length scale will increase without any bound. 
If a single growing length scale characterizes the evolving system then it is said to 
scale dynamically. Bray and Rutenberg [3j proposed a very general technique, known 
as the energy scaling approach, to estimate growth laws in purely dissipating systems 
that scale dynamically. However their scheme could also be applied to find out the 
relation between various length scales for a system in which dynamical scaling does 
not hold. If the growth law observed is different from their estimation then we can 
say that the system violates dynamical scaling. There are a large number of systems 
where the failure of dynamical scaling is observed, e.g. one-dimensional XY model 
|3], the nonconserved two-dimensional 0(3) model [H], the conserved spherical model 
etc. 

The model we have studied is described by the Hamiltonian, 

# = - £ {<t>i,<t>i? 

<i,j> 

where (ft is the usual 0(n) vector spins. Due to the spin inversion symmetry, the 
model represents a uniaxial nematic. The phase ordering of the same model was 
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studied by Blundell and Bray using a Cell Dynamic Scheme for d=2, n=2 and 
d=3, n=3 and in the present work we have studied it for d=2 and n=3. It differs 
from the usual two-dimensional 0(3) model due to its local inversion symmetry. 
While in the two-dimensional 0(3) model there is no stable topological singular 
defect, in the present model, due to presence of local inversion symmetry, the order 
parameter space, instead of being a simple three dimensional sphere (as in case of 
usual 0(3) models), is a three dimensional sphere with antipodal points identified. 
This gives rise to topologically stable singular point defects of strength ±1/2 (where 
the director rotates around the defect core by 180°). The mapping of other 1/2 
integrals defects are homotopically equivalent to the mapping of 1/2-defects. A 
— 1/2 defect configuration is continuously deformable to a 1/2 defect configuration. 
Integral defects are topologically unstable due to the so called 'escape to the third 
dimension'. Hence the first or the fundamental homotopy group of the system is just 
the two element group Z 2 ({0,1}) [B|. 

In general, the 0(n) model with d=n-l supports nonsingular topologically stable 
extended spin configurations carrying integral topological charge, known as topo- 
logical textures (or anti-textures, for negative topological charge )[SJ. In the two- 
dimensional 0(3) model the textures are known as skyrmions, instantons or baby 
skyrmions. The various length scales associated with these weakly interacting tex- 
tures evolve with different growth laws in the one-dimensional XY model and in the 
two-dimensional 0(3) model, which give rise to scaling violation in these systems 
[HE]. In the two-dimensional 0(3) model the minimum energy configuration for 
an isolated texture is obtained by stereographically projecting the order parameter 
sphere on the physical space [Hj. The configuration covers the order parameter space 
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exactly once and hence the texture is associated with a topological charge 1. In 
the present model one would expect the presence of the two-dimensional 0(3) like 
textures, but our effort to find out them using the algorithm prescribed by Berg and 
Luscher [lOjresulted in detecting no textures at all. This may be explained purely 
on the basis of homology of the order parameter space. Hind marsh [XT] on the ba- 
sis of topology or more specifically homology of the order parameter space (which 
is the projective plane RP 2 in the nematics) have shown that in three dimensional 
quenched nematics the probability of occurrence of monopoles is very low. Unlike in 
the Heisenberg model, in real nematics in order to get a monopole the order parame- 
ter space has to be covered twice and a special arrangement over many uncorrelated 
domains is required. This is responsible for a very low probability (~10 -8 ) of occur- 
rence of the monopoles. In case of the two-dimensional RP 2 model (an example of 
which is the present model) similar argument should also be valid for the textures 
and this perhaps explains why we could not find the textures in this model. In case of 
the two-dimensional 0(3) model the different growth rates associated with internal 
and external length scales of the extended textures are responsible for the failure of 
single length scaling [El El- Since in the present model textures (or antitextures) are 
highly suppressed due to topological reason, the scaling violation is less likely. In the 
present paper we have established that the system scales dynamically. 

We have used Cell Dynamic Scheme[3 H3] for studying coarsening dynamics of 
the soft spin version of the concerned model. The discrete time updating relation is, 

+ E (f) n (i) tanh(\(f) n (i)\) 



(f>n+i(i) = D 



J2(<f>n(i), <f>n(j)) <f>n(j) ~ <Pn{i) 
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Figure 1: Scaling plot of the correlation function C(r,t) against r/L cor (t) for a 256x256 
lattice (D=0.1, E=l.l) obtained after collapsing the correlation function at different 
time steps (as indicated in figure3). The correlation length is obtained by using C 
(L cor (t),t) = 0.3. The agreement of Bray Puri prediction [19] for 0(2) (• • ■) model 
with the scaled correlation function(+) (t=400) is shown in the inset. The BP 
function f BP (x) = (1/tt) exp(-x 2 /2) [B(l/2, 3/2) f F(l/2, 1/2, 2; exp(-x 2 )). The 
maximum value of the L cor obtained is 19.35 which is much smaller the linear size of 
the lattice, i.e. 256. 

The sum is over nearest neighbors of i. The parameter D is called the diffusion con- 
stant, which determines the strength of the coupling between various cells evolving 
with time. The value of the parameter E should always be greater that unity and 
it determines the depth of quench [T3]. In the above discrete time updating relation 
the unit vectors (represented by hats) are used for stability of the iteration process. 
However one must avoid using both 0(j) as unit vectors as this leads to a freezing of 
the configuration in some metastable region [Zj. 

The phase ordering kinetics of the two-dimensional uniaxial nematic have been 
studied in details by Zapotocky et al. [T2|. Using a Cell Dynamic Scheme, they have 
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Figure 2: Log-log plot of the structure factor scaling function of a 256x256 lattice 
(D=0.1, E=l.l). The first moment <k> is used for rescaling of momenta. The 
straight line (dotted) has a slope of -4, which indicates the validity of the generalized 
Porod's Law. 

shown that dynamical scaling is violated in two-dimensional uniaxial nematics films. 
They observed different values of the growth exponents, in the familiar algebraic 
growth law L(t) ~ (0 is known as growth exponent), corresponding to different 
length scales. In determining the effective growth exponents they used the time 
range between 200 and 2000. However, as indicated by Rojas and Rutenberg [16J, 
in context of the issue of dynamical scaling in two-dimensional XY model quenched 
from above T KT (the Kosterlitz - Thouless transition temperature), that in order to 
decide whether a system violates dynamical scaling or not, one must find the effective 
growth exponent in the true asymptotic limit after it is constant with time and before 
the finite size effect starts playing its role. They observed no violation in dynamical 
scaling in the two-dimensional XY model. Like integral singular point defects (known 
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Figure 3: The attempted collapse of the correlation function for a 256x256 lattice 

1 /2 

(D=0.1 and E=l.l) with respect to the the defect separation length La e f = l/p de f. 
The failure of collapse at the initial stages of phase ordering does not indicate the 
violation of dynamical scaling. In the asymptotic limit the correlation function scales 
well with respect to the defect separation, indicating the proportionality of L cor (t) 
and L def (t). 
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Figure 4: The plot of t/L 2 (t) vs. ln(t) for three lengths L cor (x), L de j (+) and L k 
(*) for a 512x512 lattice (D=0.1 and E=l.l). The observed linear dependence at 
late times (over a wide range from t=5500 to t=14400) indicates that the dynamical 
scaling growth law L(t) ~ (t//n(t/t )) 1,/2 holds. However the time scale t is found 
to be non-universal. 
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as vortices) present in the two-dimensional XY model, the present two-dimensional 
model also supports topologically stable ±1/2 disclination points. So it is expected 
that the growth laws should be similar and this is the main finding of the present 
work. In figure 5, we have shown the Schlieren patterns in a 180x180 uniaxial nematic 
placed between crossed polarizers at different times indicated in the figure. The 
patterns were obtained in the same way as discussed in references [121 EH EH] ■ By 
performing the simulation for sufficiently long run to get the true asymptotic limit, 
in the present work we have shown that, the two-dimensional uniaxial nematic scales 
dynamically by establishing that the same asymptotic growth law is valid for various 
length scales. Instead of the usual growth law, the system was found to scale 
asymptotically in a manner similar to the two-dimensional XY model quenched from 
above Tkt with the growth law, L(t) ~ (t/ln(t/t )Y^ 2 ( £ nonuniversal time scale 
) [101 Ej. We have performed our simulation with two lattices sizes 256x256 and 
512x512. By comparing the results of the two lattice sizes, we did not find any 
significant finite size effect upto the time limit we have investigated. 
The normalized correlation function in the present model is given by, 

C(r,t) = 3/2<(0(O),0(r)) 2 >-l/2 

where < > represents the average over various random initial states (random length 
and magnitude). The scaling form of correlation function is given by, 

C(r,t) = f(r/L cor (t)) 
Where the L cor (t) is the length scale required to collapse the correlation functions 
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for different time. In the Figurel we have shown the scaling plot of C(r,t) averaged 

over 20 initial states for a 256x256 lattice. 

The structure factor scales with respect to L\. — 1/ < k >|16|. where <k> = 
S(k,t) Kj J2S{k,t) , is the first moment of structure factor. The scaling form of 

the the structure factor is given by, 

S(k,t) = L 2 k g(kL k (t)) 

In Figure2 we have shown the plot of ln(Lfc~ 2 S(k,t)) against (kL^). From general- 
ized Porod's law the large-momentum structure factor for a two-dimensional system 
with point defects should be proportional to pd e j k~ A [TJ [121 HU [THJ. In the present 
system the density of the point defects (pdef) scales as L^/ -2 , where Ld e f is the 
typical defect separation length. In the large-momentum limit we obtained the the 
slope of ln(Lfc~ 2 S(k,t)) versus kLfc plot equal to —4 as shown in the Figure2, which 
verifies Porod's Law. The good collapse of the tail, verifies that and La e f have 
asymptotically the same growth law. In Figure3, where we have tried the collapse of 
correlation function with respect to the defect separation, poor collapse at the initial 
stage of dynamics however does not indicate the violation of dynamical scaling. At 
late stages the collapse shown in Figure3 is quite acceptable. In Figure4 we have 
verified that the growth law for all concerned lengths are asymptotically same as 
that for the two-dimensional XY model, because of the linear dependence of t/L 2 (t) 
on ln(t) . Figure4 also verifies that the present system does not violate dynamical 
scaling. The correlation function that we obtained in our simulation also tallies with 
the Bray Puri prediction |19J of the equal time correlation function for 0(2) model, 



10 



i.e. 



C(r,t) = (jM[B(1/2, 3/2)] 2 F(l/2, 1/2 , 2 ; 7 2 ) 

where, 

j(r,t) = exp(—r 2 /8t) 

and B(x,y) is the beta function, and F(a, b, c; z) is the hypergeometric function. 

Thus in the scaling form we have C(r,t)=fap(x) (BP stands for Bray and Puri), 
where x=r/L(t) and L(t) = (4t) 1/ ' 2 . The logarithmic factor is not correctly reflected 
in that function. However for comparison we have plotted the fap(x) and the scaled 
correlation function (for t=400) in the inset of Figurel. We have performed our 
simulation with higher values of D(0.5) and got similar asymptotic results. Higher 
values of D are useful in achieving the asymptotic regime faster. However the finite 
size effect is also more prominent in case of large D. We have not considered noise 
in the time evolution equation, hence are effectively working at T=0. However it 
is known that quenching to T=0 may lead to metastable freezing [20]. In order to 
check that our results are not influenced by such freezing, we performed a number 
of simulations (almost 100 steps) with noise, by adding a constant amplitude (of 
the order of 0.1) random configuration to the order parameter. The noise amplitude 
used was enough for generating large number of pairs of disclination point defects. 
We could not find any discrepancy with the results obtained without noise. 

To summarize we would like to focus on the main findings of the paper. We 
have confirmed that dynamical scaling is not violated in two-dimensional nematic 
with order parameter dimensionality three and asymptotically the growth laws are 
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same as that of the two-dimensional XY model quenched from above Tkt (i.e. the 
initial state with free vortices). Consideration of topological defect in the issue of 
dynamics is very necessary, because the structure of the defects determines the large- 
momentum dependence of the structure factor, which has an important role in the 
determination of the growth laws (3). In the O(n) model with n<2, the topological 
defects dominate the dynamics. Since both the present two-dimensional model and 
the two-dimensional XY model support singular point defects, it is expected that 
the dynamics should be similar and this is established in the present work. In the 
present work we were able to establish the expected result and achieved the true 
asymptotic limit of the dynamical scaling by performing longer simulation with the 
help of increased computational power now available. 
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